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Abstract. We study numerically the influence of strong Coulomb repulsion on the current through molec- 
ular wires that are driven by external electromagnetic fields. The molecule is described by a tight-binding 
model whose first and last site is coupled to a respective lead. The leads are eliminated within a pertur- 
bation theory yielding a master equation for the wire. The decomposition into a Floquet basis enables an 
efficient treatment of the driving field. For the electronic excitations in bridged molecular wires, we find 
that strong Coulomb repulsion significantly sharpens resonance peaks which broaden again with increasing 
temperature. By contrast, it has only a small influence on effects like non-adiabatic electron pumping and 
coherent current suppression. 

PACS. 05.60.Gg Quantum transport - 85.65.-|-h Molecular electronic devices - 72.40.-|-w Photoconduction 
and photovoltaic effects - 73.63.-b Electronic transport in mesoscopic or nanoscale materials and structures 



1 Introduction 

Recent experiments on the conductance of single organic 
molecules opened a new direction in mesoscopic trans- 
port [1-3] . Of particular interest is thereby the influence of 
electronic and vibronic excitations of the molecules which 
leave their fingerprints in the resulting current-voltage 
characteristics. At low temperatures, these effects become 
more pronounced [4] . Much of our knowledge about excita- 
tions of molecules is based on spectroscopy, i.e. the optical 
response to light. In the context of molecular conduction, 
it has been proposed to study as well the signatures of 
such excitations in the transport quantifiers like the cur- 
rent [5-7] and its fluctuations [8] . Such experiments are at 
present attempted, but clearcut evidence for the proposed 
effects is still missing because the irradiation also causes 
unwanted thermal effects in the contacts, which in today's 
setups seem to dominate. One possibility to protect the 
contacts against the radiation is using the evanescent light 
at the tip of a near-field optical microscope. 

Coupled quantum dots represent a setup with proper- 
ties similar to those of molecular wires [9,10], albeit at dif- 
ferent length and energy scales. As compared to molecular 
wires, they are more stable and tunable, but have the dis- 
advantage that only a few dots can be coupled coherently. 
The transport properties of these "artificial molecules" 
can be significantly modified by microwaves [10]. It has 
for example been demonstrated experimentally that res- 
onant excitations between the levels of double quantum 
dots result in the so-called photon-assisted transport, i.e. 
a significant enhancement of the dc current [11, 12]. 

In addition to photon-assisted transport [13], other less 
intuitive phenomena have been predicted in this context. 
For example the so-called coherent suppression of tun- 



nelling in a double- well potential due to a time-dependent 
dipole force can also be found in ac driven transport through 
coupled quantum dots: For characteristic ratios between 
the amplitude and the frequency of the driving, it has 
been predicted that the dc current [14-16] and the shot 
noise [17] will be suppressed. A further prominent effect 
is adiabatic electron pumping, which is the generation of 
a dc current by means of a periodic variation of the con- 
ductor parameters in the absence of any net bias [18-21]. 
It has been proposed [22,23] and experimentally demon- 
strated [10,11] that pumping is more effective at internal 
resonances, i.e., beyond the adiabatic limit where, in addi- 
tion, the pump current possesses a surprisingly low noise 
level [24]. 

Periodically time-dependent quantum systems can be 
described very efRciently within a Floquet theory which 
originally has been derived for driven closed quantum sys- 
tems [25] and later been generalised to dissipative quan- 
tum systems [26,27]. Furthermore, it is possible to derive 
Floquet theories for the description of transport through 
mesoscopic conductors which are connected to external 
leads. For cases in which electron-electron interactions 
do not play any role, one can derive a Floquet scatter- 
ing theory that provides exact expressions for the cur- 
rent [17,28] and its noise [8,17]. Treating the coupling 
between the conductor and the leads perturbatively, one 
can obtain a master equation for the reduced density oper- 
ator of the wire. This enables a rather efficient treatment 
of time-dependent transport after decomposing the wire 
density operator into a Floquet basis. Then it is possi- 
ble to study relatively large driven conductors [29] and 
to include also electron-electron [30] and electron-phonon 
interactions [31]. If the time-dependent field consists of 
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Fig. 1. Bridged molecular wire model consisting of iV = 5 
sites with internal tunnelling matrix elements A and effective 
wire-lead coupling strengths /l/r- 

one or a few laser pulses, it is possible to obtain the den- 
sity operator by propagating the Liouville-von Neumann 
equation [32]. 

In this work, wo derive in Sections 2 and 3 a Flo- 
quet master equation formalism that captures situations 
in which strong Coulomb repulsion restricts the excess 
charge residing on the conductor to a single electron. For 
later reference, we adapt in Section 4 our approach to 
the case of spinless electrons and non-interacting elec- 
trons. We then investigate in Section 5 the role of strong 
Coulomb repulsion for photon-assisted transport through 
bridged molecular wires, non-adiabatic electron pumping, 
and coherent current suppression. 



2 The wire-lead model 

The system of the driven wire, the leads, and the coupling 
between the molecule and the leads, as sketched in Fig. 1, 
is described by the Hamiltonian 



Ti-it) — Hwiro(i) + Wlcads + Wwiro-loads- 



(1) 



The wire is modelled within a tight-binding description by 
the molecular orbitals |n), n = 1, . . . , A'', so that 

^wire(^) ~ ^ ^ ^nn'i^^ ^^g^Yifs~^ ~^^wiTe{-^wiTe 1)5 (2) 

where cj^g (c„s) creates (annihilates) an electron with spin 
s at site |n) and [c„,,, cj^,j,,]+ = (5n.ri''5s,s'- The influence 
of a driving field entails a T-periodic time-dependence on 
the single-particle Hamiltonian iJ„„/(i) = Hnn'{t + T). 
The last term in Eq. (2) captures the electron-electron 
interaction within a capacitor model and the operator 
■^wire = s ^ns'^ns describes the number of excess elec- 
trons residing on the molecule. Below we shall assume that 
U is so large that only states with zero or one excess elec- 
tron play a role. 

The first and the last site of the molecule, |1) and \N), 
couple via the tunnelling Hamiltonian 

Wwire-leads = ^{Vhq clg^C^^, + Fr, ^ggCjvJ + H-C (3) 
q,s 



to the respective lead. The operator Cl^^, (cr^s) creates an 
electron in the left (right) lead in the state \ljqs) which is 
orthogonal to all wire states. It will turn out that the in- 
fluence of the tunnelling Hamiltonian is fully characterised 
by the spectral density 



re{€) = 27rJ2\Ve<ifS{€ 



(4) 



which becomes a continuous function of the energy e if 
the lead states are dense. If all relevant lead states are 
located in the centre of the conduction band, the energy- 
dependence of the spectral densities is not relevant so that 
they can be replaced by a constant, /^(e) = F^. This de- 
fines the so-called wide-band limit. 

The leads are modelled as ideal Fermi gases 



■^leads — ^ (Cg <^qs^hqs + ^q ^Rqs^Kqs) ' 



(5) 



q,s 



which are initially at thermal equilibrium with the chemi- 
cal potential /Ul/r and, thus, are described by the density 
operator 

Pleads.eq OC Cxp [-(Wfeads - ML-^Vl - IJ.nNn)/kBT], (6) 

where = J2q c^pCqi denotes the electron number in lead 
i = L, R. Then all lead properties can be expressed in 
terms of the expectation value 



{^qs^^t'q's') = Su'Sqq'Sss'MCq), 



(7) 



where fi{e) = (I +e^^-^'''>/''^'^)-'^ denotes the Fermi func- 
tion. Since a typical metal screens all electric fields with 
a frequency below the plasma frequency, we assume that 
the bulk properties of the leads are not affected by the 
laser irradiation. 



3 Master equation approach for strong 
Coulomb repulsion 

3.1 Perturbation theory 

Most master equation approaches to electron transport 
[30 34] are based on a perturbative treatment of the wire- 
lead Hamiltonian i?wire-ieads- Within these approaches, it 
is possible to include within a Floquet theory the response 
to a time-periodic field exactly [8,35]. The derivation of 
a master equation starts from the Liouville-von Neumann 
equation ihp{t) = [H{t),p{t)] for the total density opera- 
tor p{t) for which one obtains by standard techniques the 
approximate equation of motion 



P{t) = - ^[Hwiie{t) + i?leads,P(t)] 



(8) 



1 r°° 



i-leads) [-^Avire- leads 

{t-T,t),p{t)]]. 
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Here the first term corresponds to the coherent dynam- 
ics of both the wire electrons and the lead electrons, while 
the second term describes resonant electron tunnelling be- 
tween the leads and the wire. The tilde denotes operators 
in the interaction picture with respect to the molecule and 
the lead Hamiltonian without the molecule-lead coupling, 
X{t,t') = U^{t,t') XUa{tJ/), where Uq is the propagator 
without the coupling. The net (incoming minus outgoing) 
electrical current through the left contact is given by mi- 
nus the time-derivative of the electron number in the left 
lead multiplied by the electron charge — e. From Eq. (8) 
follows for the current in the wide-band limit the expres- 
sion 

X {{c\c,{t, t - r))/L(e) - {c,{t, t - T)c\)Me)}, 

where fe = i — fi and (• • • ) = trwirc Pwirc • • • denotes the 
average with respect to the wire density operator, which 
still has to be determined. We emphasise that the expec- 
tation values in Eq. (9) depend directly on the dynamics 
of the isolated wire and are thus influenced by the driving. 

3.2 Floquet theory 

An important feature of the current formula (9) is its de- 
pendence on solely the wire operators while all lead opera- 
tors have been eliminated. Therefore it is sufHcient to con- 
sider the reduced density operator pwire = trieads P of the 
wire electrons. However, p^ic still obeys a time-dependent 
master equation whose direct solution requires quite some 
effort, in particular if one is interested in the behaviour 
at long times. This effort can be reduced significantly by 
exploiting the fact that the master equation (8) inher- 
ited from the total Hamiltonian Ti.{t) a T-pcriodic time- 
dependence, which opens the way for a Floquet treatment. 

3.2.1 Fermionic Floquet operators 

One possibility for such a treatment is to use the Floquet 

states of the central system, i.e. the driven wire, as a ba- 
sis. Thereby we also use the fact that in the wire Hamilto- 
nian (2), the single-particle contribution commutes with 
the interaction term and, thus, these two Hamiltonians 
possess a complete set of common eigenstates. Here we 
start by diagonalising the first part of the Hamiltonian 
which describes the single-particle dynamics determined 
by the time-periodic matrix elements i?„„'(t). According 
to the Floquet theorem, the corresponding (single parti- 
cle) Schrodinger equation possesses a complete solution of 
the form 

|'^'o:(t))=e-°*/''|^„(t)), (10) 

with the so-called quasienergies and the T-periodic Flo- 
quet states 

fe 



The Floquet states and eigenenergies are obtained by solv- 
ing the eigenvalue problem 

(^^\n)Hr,n'{t){n'\-ih^)\Mt))=eaM)), (12) 

n,n' 

whose solution allows one to construct via Slater deter- 
minants many-particle Floquet states. In analogy to the 
quasimomcnta in Bloch theory for spatially periodic po- 
tentials, the quasienergies come in classes 

Ea.fe = + khO, k gZ, (13) 

of which all members represent the same physical solution 
of the Schrodinger equation. Thus we can restrict ourselves 
to states within one Brillouin zone like for example < 

Cq < h[2. 

For the computation of the current it is convenient to 
have an explicit expression for the interaction picture rep- 
resentation of the wire operators. It can be obtained from 
the (fermionic) Floquet creation and annihilation opera- 
tors [36] defined via the transformation 

Casit) =^{(Pa{t)\n)Cns- (14) 
n 

The inverse transformation 

Cns = ^{n\<fia{t))Cas{t) (15) 
a 

follows from the mutual orthogonality and the complete- 
ness of the Floquet states at equal times [27]. Note that 
the right-hand side of Eq. (15) becomes time independent 
after the summation. 

The Floquet annihilation operator (14) has the inter- 
action picture representation 

Cas{t,t') = Ul{t,t/)Cas{t)Uo{t,t') (16) 
^ g-i(e„-FC/iV^.„)(t-t')/ftc^^(i')^ (17) 

with the important feature that the time difference t — t' 

enters only via the exponential prefactor. This will allow 
us to evaluate the r-integration of the master equation (8) 
after a Floquet decomposition. Relation (17) can easily be 
shown by computing the time derivative with respect to t 
which by use of the Floquet equation (12) becomes 

^Cas{t,t') = -U<^a + UN^ire)Cc,s{t,t'). (18) 

at n 

Together with the initial condition Ca{t',t') = Ca{t') fol- 
lows relation (17). Note that the time evolution induced 
by Hwire(0 conserves the number of electrons on the wire. 

3.2.2 Master equation and current formula 

In order to make use of the Floquet ansatz, we decompose 
the master equation (8) and the current formula (9) into 
the Floquet basis derived in the last subsection. For that 
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purpose we use the fact that wc are finally interested in 
the current at asymptotically large times in the limit of a 
large interaction U. The latter has the consequence that 
only wire states with at most one excess electron play a 
role, so that the density operator Pwire can be decomposed 
into the 2 A'' + 1 dimensional basis {|0), c1jj,(t) |0)}, where 
|0) denotes the wire state in the absence of excess electrons 
and s =t, J.- Moreover, it can be shown that at large times, 
the density operator of the wire becomes diagonal in the 
electron number A^wire- Therefore a proper ansatz reads 



Pwire(i) = |0)P00(i)(0| + ci,\O)pas,Ps'{t){O\C0s'- 

(19) 

Note that we keep terms with a ^ (3, which means that we 
work beyond a rotating-wave approximation. Indeed in a 
non-equilibrium situation, the off-diagonal density matrix 
elements pa/s will not vanish and neglecting them might 
lead to artefacts [8,37]. In the context of molecular wires, 
such a treatment of strong Coulomb repulsion by a re- 
striction to at most one excess electron has recently also 
been applied to incoherent [38^0] as well as to coherent 
transport [41]. 

By inserting the decomposition (19) into the master 
equation (8), we obtain an equation of motion for the ma- 
trix elements Pasj3s' — (OlcasPwireC^g/ 10). We evaluate the 
trace over the lead states and compute the matrix element 
(0|cas(i) ■ ■ ■ cjjg, (i)|0). Thereby we neglect the two- particle 

terms which are of the structure c^j^c^jglO) (0|c/3sCq,s. For- 
mally, these terms drop out in the limit of strong Cou- 
lomb repulsion because they are accompanied by a rapidly 
oscillating phase factor ex]){—iUN^iieT/fi). Then the r- 
integration results in a factor f-L{ea,k + U) which vanishes 
in the limit of large U. Since the total Hamiltonian (1) is 
diagonal in the spin index s, we find that the density ma- 
trix elements pas, 0s' are spin-independent as well so that 
after a transient stage 



PaT,/3T(0 = PaUlit) = Pap{t) (20) 



and Pa^.fii = 0. Moreover, at large times, the density op- 
erator (19) will acquire the time periodicity of the driving 
field [8] and, thus, can be decomposed into the Fourier 
series 

Pa/3(i) = E^"''''W,fe (21) 
k 



and pQo{t) accordingly. 

After some algebra, we arrive at a set of N"^ coupled 
equations of motion for pai3{t) which in Fourier represen- 



tation read 
i(ea - e/3 - khfi)pai3,k 

= Yj ('^a>'='+fe"|l)(l|V/3,fc-Ffc")P00,fe' 
k',k" 

X (/l( ) + fhief3,k+k")) 

2~ Yj {'Pa,k'+k"\^){MVa',k+k>')pa>l3,k> fl,{£a',k+k") 

a',k',k" 

- {.'^l3',k'+k"\'^)i}-W0,k+k")Pal}',k' h{^0',k'+k") 

I3',k',k" 

+ same terms with the replacement 1, L ^ A'', R. 

(22) 

In order to solve these equations, we have to eliminate 
PoOj/c which is most conveniently done by inserting the 
Fourier representation of the normalisation condition 

tr Pwire {t) = Poo (i) + 2 X] Poia (*) = !• (23) 
a 

In order to obtain for the current an expression that is 
consistent with the restriction to one excess electron, we 
compute the expectation values in the current formula (9) 
with the reduced density operator (19) and insert the Flo- 
quet representation (15) of the wire operators. Performing 
an average over the driving period, we obtain for the dc 
current the expression 

,fe'-|-fc|l)(l|¥'a,fc)PQ/3,fc'/L(ea,fe) 

a, A: 0,k' 

— X^(</'a,/£'-|-fe|l)(l|¥'a,fe)pOO,fe'/L(ea,fe)) • 
k' 

(24) 

Physically, the second contribution of the current formula 
(24) describes the tunnelling of an electron from the left 
lead to the wire and, thus, is proportional to poo/l which 
denotes the probability that a lead state is occupied while 
the wire is empty. The first terms corresponds to the re- 
versed process namely the tunnelling on an electron from 
site |1) to the left lead. 

The results of this section allow us the numerical com- 
putation of the dc current through a driven conductor 
in the the following way; First, we solve the quasienergy 
equation (12) which provides the coeSicients {'^a,k \n) . Next, 
we solve the master equation (22) and insert the solution 
into the current formula (24). 

4 Separating interaction and spin 

In order to determine the role of a strong interaction, 
wc shall compare below our resiilts to the non-interacting 
case. Moreover, a particular consequence of strong Cou- 
lomb repulsion is the mutual blocking of different spin 
channels. This motivates us to also compare to the case of 
spinless electrons which is physically realised by spin po- 
larisation. In this section, we adapt our master equation 
approach to these situations. 



Franz J. Kaiser et al.: Coulomb repulsion effects in driven electron transport 



5 



4.1 Spinless electrons 

In order to describe spinless electrons, we drop in the ini- 
tial Hamiltonian all spin indices. Physically, this limit is 
realised by a sufficiently strong magnetic field that po- 
larises all electrons contributing to the transport. By the 
same calculation as in Section 3, we then obtain for the 
current also the expression (24) but without the prefactor 
2. The factor 2 is also no longer present in the normalisa- 
tion condition (23) which now reads 

tr Pwire(i) = Poo(i) + X] Pca{t) = 1- (25) 

a 

4.2 Non-interacting electrons 

In the absence of interactions, U = i), each spin degree of 
freedom can be treated separately. Still one has to consider 
for each spin projection up to A'' electrons so that the 
relevant Hilbcrt space has the dimension 2^. Therefore, 
it is more efficient to consider the single-particle density 
matrix 

. (26) 

which is of dimension N'^ and nevertheless contains all 
relevant information. The Fourier decomposition in the 
last expression uses the fact that, at asymptotically large 
times, Rap{t) becomes time-periodic. 

We express the time derivative of Rapit) with the mas- 
ter equation (8) and insert for Ca{t,t — r) the Floquet 
representation (18). After some algebra, we obtain for the 
Fourier coefficients Ra0,k the equation 

i(ea - e/3 + khf2)Rai3,k 

= -y E ( E {'^l3,k'+k"\i){'i-\f0',k+k")Rafl',k' 
k' P',k" 

+ E {'^oc' ,k'+k"\'^){'^\'^a,k+k")Ra'0,k' 

a',k" (27) 

-{'Pl3,k'-k\i){M'Pa,k'}fh{ea,k') 
-(<P/3,fe'|l)(l|Va,fe'-|-fe)/L(e/3,fe')) 

-|- terms with the replacement L, 1 ^ R, N. 

In contrast to the master equation (22) for the limit of 
strong Coulomb repulsion, no blocking factors 1 — /l/r 
emerge. This is characteristic for the non-interacting limit 
and is also found in the exact scattering formula [8]. 

Using the single-particle density matrix Rap, one can 
evaluate the current expectation value (9). Thereby the 
spin degree of freedom enters simply as a factor 2 so that 
we obtain for the current per spin the expression 

1= [^y^AV0,k'+k\'i-)('i-\^a,k)Ra0,k' 

a,k 0,k' (28) 

- |(lbc«,fc)|VL(ea,fc) • 



For a detailed derivation see Ref. [36]. 

5 Interplay of dipole radiation and Coulomb 
repulsion 

In our model Hamiltonian (1) we have already specified 
the interaction, the lead Hamiltonian, and the wire-lead 
coupling. By contrast, for the Hamiltonian of the driven 
wire, we have thus far only assumed that the external field 
is periodic in time. In the following we focus on models 
where the single-particle dynamics is determined by the 
N-site tight-binding Hamiltonian 

N-l 

Hnn'{t) = - ^ E + 1| + I" + 1>("I) 

(29) 

+ ^{En + Axn cos{Qt)}\n){n\. 

n=l 

Neighbouring sites are coupled by a tunnel matrix element 
A. The onsite energies i?„ are modulated by a harmon- 
ically time-dependent dipole force, where the amplitude 
A is given by the electrical field amplitude multiplied by 
the electron charge and the distance between neighbour- 
ing sites. Xn = \{N +l — 2n) denotes the scaled position of 
site \n). Depending on the onsite energies En. one observes 
various phenomena which we discuss in the following. 

5.1 Resonant excitations of bridged molecular wires 

A frequently studied model is the so-called bridged molec- 
ular wire [42] sketched in Fig. 1. It consists of a molecule 
with N sites, where the first and the last site — the donor 
|1) and the acceptor \N) — are connected to respective 
leads. The energies of these two sites are assumed to be 
close to the chemical potentials of the respective leads, 
Hi^ < El = En < Mr- The remaining A''— 2 orbitals lie well 
above the chemical potentials at an energy Eb ^ A which 
defines the bridge height. In the presence of laser excita- 
tions, we expect an enhanced current whenever the energy 
quanta of the laser field HH matches the energy difference 
between donor/acceptor and one of the N—2 bridge levels. 
This photon-assisted tunnelling indeed exhibits resonance 
peaks which obey the scaling law /peak oc A"^ /{N — l)r [8]. 

Figure 2 shows the dc current as a function of the laser 
frequency for the three approaches considered herein. In 
all these cases, the resonance peaks are at the same fre- 
quencies. As a further interesting feature, we find that in 
the strongly interacting case U = oo, the current does not 
depend on whether one includes the spin degree of free- 
dom. This behaviour can be understood qualitatively in 
the following way: If there is no excess electron on the wire, 
an electron can enter from both the "spin- up channel" and 
the "spin-down channel". By contrast, when leaving the 
wire, the electron spin is preserved so that only either 
channel is available. This means that including the elec- 
tron spin effectively modifies the in-tunnelling rates by 
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Fig. 2. Photon-assisted current in a bridged molecular wire 
with length N = 8 and height Eb = WA. Leads with a chemi- 
cal potential difference eV = 5A couple to the wire with am ef- 
fective strength 7^ = 0.1^. The driving amplitude is j4 = 0.3A. 
For U — oo, the current for spinless electrons coincides with 
the one for real electrons. 



a factor 2. The electron dynamics within the wire, how- 
ever, is described by the coherent first term of the master 
equation (8) which is spin-independent. Consequently, we 
expect that the spin only plays a minor role whenever the 
tunnelling from the donor to the acceptor becomes the 
bottleneck for the electrons. This is indeed the case for 
the transport across a barrier considered here. 

As compared to the non-interacting case, one notices 
that strong Coulomb repulsion modifies the shape of the 
peaks: They become slightly higher and much sharper. 
This effect is more pronounced for large wires. A closer 
look at the resonance peaks for wires with up to = 10 
sites (not shown) indicates that the resonance widths scale 
roughly like 1/N. This is possibly caused by the fact that 
the Coulomb repulsion reduces the number of available 
(many-particle) wire states and, thus, the number of decay 
channels. This relates to the observation made in Ref. [41], 
namely that Coulomb repulsion can improve quantum co- 
herence and thereby enhance the current. 

Since quantum coherence is also temperature depen- 
dent, it is natural to ask whether the resonance peaks 
become sharper for lower temperatures. Figure 3 shows 
the width of the central peak of Fig. 2 as a fimction of the 
temperature. While in the non-interacting case, the peak 
widths are essentially temperature independent, the sit- 
uation changes for strong Coulomb repulsion: There one 
finds that with an increasing temperature, the peaks be- 
come roughly twice as broad once the temperature exceeds 
T = A/k-Q. We attribute this behaviour to the reduced 
coherence for thermal energies that are larger than the 
tunnelling matrix element. 



5.2 Non-adiabatic electron pumping 

Another well studied phenomenon in driven transport is 
coherent electron pumping, i.e., the creation of a non- 
vanishing dc current by ac fields in the absence of any net 
bias. For adiabatically slow driving, this effect exists only 
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Fig. 3. Temperature dependence of the widths of the reso- 
nance peaks at Q ~ Q.TlA/h. The graph shows the full widths 
at half of the maxinniin. All other parameters arc as in Fig. 2. 
For U = oo, the values for spinless electrons and real electrons 
coincide. 
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|2) 



|1> 



MR 



Fig. 4. Tight-binding model for two coupled quantum dots 
in pump configuration, i.e. in the absence of a bias voltage 
but with an internal bias E2 — Ei ^ which breaks reflection 
symmetry. 



in the absence of time-reversal symmetry [18-20]. Beyond 
the adiabatic regime, this is no longer the case: For fast, 
time-periodic driving fields, it can be shown that the rele- 
vant symmetry is the so-called generalized parity which is 
defined as the invariance under spatial reflection in com- 
bination with a time shift by half a driving period [8]. 
Non-adiabatic electron pumping is particularly interesting 
because at internal resonances of the central system the 
pump current can assume rather large values [10,22,35], 
while at the same time the current noise is remarkably 
low [24]. 

For studying the influence of strong Coulomb repulsion 
on non-adiabatic electron pumping, we consider the setup 
sketched in Fig. 4. Of particular interest is the parameter 
regime with large internal bias and intermediate dot-lead 
coupling because in this regime, the current-to-noise ratio 
is most favourable [24] . The currents obtained for the three 
considered approaches are shown in Fig. 5. Again we find 
that the spin degree of freedom is not of major influence, 
which indicates that the transport is governed by internal 
excitations; cf. the discussion in the preceeding subsection. 

The influence of the Coulomb repulsion is a modifica- 
tion of the current peak height up to 5%. This means that 
interactions are here much less important than for photon- 
assisted transport: The reason for this is that for our pump 
configuration, one energy level lies below the Fermi energy 
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Fig. 5. Pump current for strong Coulomb repulsion and for 
non-interacting electrons as a function of the frequency. The 
dot levels with energies Ei = —2.5A and E2 = 2.5A are couple 
to the leads with strength F = 0.3A. The driving amplitude is 
A = 3.7 A and the temperature ksT = O.OOSZi. 



while the other lies well above. Thus in equilibrium for a 
sufficiently small dot-lead coupling, the left site is occu- 
pied while the right site is empty, whatever the interaction 
strength. Thus, the double dot is populated with only one 
electron so that interactions become irrelevant. Unless the 
driving amplitude is huge, this occupation is altered only 
slightly. Consequently interactions do not modify the cur- 
rent significantly. We emphasise that for strong dot-lead 
coupling r and finite interaction U, these arguments no 
longer hold true. 



5.3 Coherent current control 

An intriguing example of quantum control is the so-called 
coherent destruction of tunnelling in a double-well po- 
tential by a suitable driving field [43], which can be ex- 
plained within a rotating-wave approximation: In the driv- 
ing dominated regime, the tunnel matrix element is es- 
sentially replaced by an effective tunnel matrix element 
Aeff = Jo{A/hS7), where Jq is the zeroth order Bessel 
function of the first kind and A the driving amplitude 
[44, 45] . Related effects have been predicted for driven 
tight-binding lattices [46]. For driven transport between 
two leads, the corresponding situation has been investi- 
gated only recently: It has been found that driving fields 
that suppress tunnelling in a closed driven system, also 
suppress the current through the corresponding open sys- 
tem [14,47]. 

A setup in which this coherent current suppression can 
be studied is sketched in Fig. 6. It is characterized by the 
fact that all internal levels lie within the voltage window, 
i.e. below one chemical potential and above the other. For 
this system, strong interaction has already a significant 
influence on the current in the absence of any driving field 
[41]: For U = 0, the system is half filled, which means it is 
populated by N/2 electron, and the current is independent 
of N. By contrast for U = 00, the stationary population 
is N/ {N + 1) and the current becomes / cx 1/(A'^ -|- 1). 



eV 



Ml 



A 



A 



|1) |2) 



|3) 



Mr 



Fig. 6. Triple quantum dot configuration for coherent current 
control: A large bias voltage and the unbiased dot levels with 
E„ = ensure optimal transport in the absence of the driving. 
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Fig. 7. Coherent current suppression as a function of the driv- 
ing amplitude A for a wire that consists of iV = 3 sites. The 

applied voltage is eV — 50A, the driving frequency fi = 5A, 
and the wire-lead coupling F = 0.5A. 



Figure 7 shows the influence of Coulomb repulsion on 
the current suppression studied in Ref. [17]. Independent 
of the interaction, one finds that the current almost van- 
ishes whenever the ratio A/hil matches a zero of the 
Bessel function Jq. If the driving amplitude is far from the 
values for which the current is suppressed, we observe the 
behaviour found for the static situation, namely that Cou- 
lomb repulsion reduces the current by a factor 1/(A'^ -|- 1) 
for spinless electrons [41]. If one considers the spin, this 
factor becomes 2/{2N + 1). In the vicinity of the current 
suppressions, by contrast, the influence of both the inclu- 
sion of the spin and the interaction is less pronounced. 
In this regime, the effective tunnel matrix element Ags is 
small, so that tunnelling along the wire happens at a low 
rate. This again indicates that whenever the transport is 
limited by the dynamics within the wire, the influence of 
interaction is rather small. 
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6 Conclusions 

We studied the influence of strong interaction on the trans- 
port properties of ac-drivcn coherent conductors. In par- 
ticular, we compared the strongly interacting case with the 
opposite extreme of non-interacting electrons. Moreover, 
we worked out the relevance of the spin degree of free- 
dom for weak wire-lead coupling. In our studies, we con- 
sidered three archetypical effects, namely photon-assisted 
tunnelling through bridged molecular wires, non-adiabatic 
electron pumping, and coherent current suppression. 

The most significant effect is found for photon- assisted 
tunnelling where Coulomb repulsion renders the resonance 
linewidths much sharper. Thus unfortunately, interactions 
might contribute to the difficulties in photon-assitcd tun- 
nelling experiments with molecular wires. By contrast. 
Coulomb repulsion is not too relevant for electron pump- 
ing in double quantum dots. For coherent current suppres- 
sion, the same holds true only for parameters for which 
the current is already significantly reduced. Outside this 
region, one finds that Coulomb repulsion reduces the cur- 
rent essentially in the same way as in the absence of driv- 
ing. 

The two extreme cases of zero and very strong interac- 
tion do not necessarily allow a simple interpolation. Thus, 
it is desireable to extend the present studies to finite values 
of the interaction strength, which requires the generalisa- 
tion of our formalism to at least a second excess electron. 
Moreover, the dc current is certainly not the only relevant 
quantity for the characterisation of the electron transport. 
Investigating the influence of Coulomb repulsion on, e.g., 
the current noise would complement the picture drawn 
above. 

We thank M. Strass and A. Nitzan for interesting discussions. 
This work has been supported by Deutsche Forschungsgernein- 
schaft through SFB 484 and SPP 1243. One of us (FJK) ac- 
knowledges funding by Bayerisches Staatsministerium fiir Wis- 
senschaft, Forschung und Kunst through Elitenetzwerk Bay- 
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